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Abstract. 

We construct a general parametrization of the baryon octet and decuplet mass 
operators including the three-body terms using the unit operator and the symmetry- 
breaking factors M'^ = diag (0,1,0) and = diag (0,0,1) in conjunction with the 
spin operators. Our parametrization has the minimal number of operators needed 
to describe all the octet and decuplet masses. Investigating the likely size of the 
three-body terms, we find that contributions of the three-body hypercharge splittings 
are comparable to those from the one- and two-body isospin splittings and that 
contributions of the three-body isospin splitting operators are very small. We prove 
that, in dynamical calculations, one must go to three loops to get the three-body 
terms. We also find that the suggested hierarchy of sizes for terms in the most general 
expression for baryon masses that involve multiple factors of M'^ and/or M" does not 
hold strictly for dynamical calculations in heavy baryon chiral perturbation theory: 
terms of a given order in a meson loop expansion may appear both with the expected 
factors of M'^ and M^, and with one factor more. 
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1. Introduction 

The study of baryon masses is one of tlie fundamentally important problems in 
nuclear and particle physics and has been of interest for many years. Progress on 
the modern theory began with the introduction of flavor SU(3) symmetry and of the 
electromagnetic mass relation derived by Coleman and Glashow [1]. It was later found in 
the nonrelativistic quark model [2l[3llll[5] that, by including only two-body interactions 
among the quarks in the baryons, there exist nine mass formulas (a.k.a. sum rules) 
connecting the eighteen charge states of the baryon octet and decuplet. One of them is 
the remarkably accurate Coleman-Glashow relation. The baryon octet and decuplet 
masses and the sum rules have been also investigated in various versions of chiral 
perturbation theory [HI [TJ [H [9l [TOl [HI [12] . As remarked elsewhere, the chiral expansion 
gives a complete parametrization of the static properties of the baryons when carried 
to high enough order, but no dynamical information unless the new couplings can be 
calculated in the underlying theory (see, for example, the analysis of the magnetic 
moment problem in [13]). 

It was Morpurgo who first constructed a general parametrization of baryon masses 
for the case of hypercharge breaking. His general parametrization method [H] , derived 
exactly from the QCD Lagrangian employing only few general properties 0, expresses the 
possible mass operators in terms of flavor- dependent terms proportional to powers of the 
strange-quark projection operator (denoted by ilf* in our work), and nonrelativistic- 
appearing products of Pauli spin operators cr. The results are completely general and 
relativistic. Morpurgo also considered the purely electromagnetic contributions to the 
baryon masses and showed that the known sum rules for mass splittings within isospin 
multiplets hold in the absence of three-body terms [15j. Dillon and Morpurgo later 
included the light quark mass terms in some of their work. For example, they briefly 
discussed the effects of md— m„ on the Coleman-Glashow and other baryon mass formula 
in [16] where the first order of the light quark mass difference — rriu is considered. 

As already mentioned, baryon masses have been studied extensively in chiral 
perturbation theory (ChPT). See, for example, [HI [71 [H [9l [TOl [Til [12] the references in 
those papers. The possible mass operators are customarily expressed in terms of traces of 
products of baryon effective fields and their derivatives using a matrix representation of 
the fields, and mesonic corrections to the initial operators are calculated in perturbation 
theory. This effective field theory approach is based on a well-defined low-momentum 
approximation to QCD with a structure determined by the general properties of QCD 
in the limit of small quark masses and momenta fTT], and gives a completely general 
description of the low-energy properties of the baryons. 

In [T8l [T9| [20] , we reconsidered the mass problem using the heavy-baryon version 
of chiral effective field theory (HBChPT) [2T]. We used a "quark" representation in 
which the baryon effective fields are labeled using three fiavor indices i, j, k and spin 

I The general properties are listed in the first footnote in [28] . 
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indices rather than the usual matrix representation [§|, and calculated the perturbative 
corrections in a meson-loop expansion. We found that this approach leads to a 
(relativistic) description of the possible mass operators in HBChPT in terms of simple 
flavor and spin operators which is identical to that developed by Morpurgo [H]. This 
equivalence is expected since HBChPT retains the low-momentum structure of QCD, 
and the effective-field results must be consistent with the general structure of operators 
and matrix elements in QCD, and give a complete description of those quantities at low 
energies where the effective field methods apply. HBChPT also provides a formalism 
in which one can make practical dynamical calculations perturbatively, an important 
consideration in our earlier work. 

In the present paper, we use Morpurgo's method [13] to construct a general 
parametrization of the baryon masses including three-body terms using the unit operator 
and the symmetry-breaking factors, M*^ = diag (0,1,0) and = diag (0,0,1), in 
conjunction with the spin operators. Our general expression includes all orders of the 
light quark mass difference rrid — mu which are not explicitly shown in the general 
parametrization for baryon masses constructed by Dillon and Morpurgo. Using a Gram 
matrix analysis, we establish a minimal set of 18 independent operators to describe 
18 baryon octet and decuplet mass states. The present work gives the overall setting 
without relying on HBChPT. 

Our objective in our earlier studies of mass splittings (and moments) was to see the 
extent to which the observed splittings (other than explicit quark mass terms) could be 
explained in terms of dynamical interactions in HBChPT and the added electromagnetic 
interactions [20l[22] . rather than just parametrized. Our approach was therefore mainly 
calculational, restricted to one- or two-loop corrections to the initial mass operators, and 
we did not include an explicit analysis of the properties behind the parametrization, or 
of the most general forms allowed. We found in that work that all the one- and two- 
body operators were generated in a meson loop expansion in HBChPT carried to two 
loops, but no three-body operators were generated. We will show here that one must 
go to three loops in HBChPT to get three-body terms. Our general demonstration of 
this and the examples of nontrivial three-body terms are also new. 

We will also investigate the likely sizes of the coefficients of the various operators, 
particularly the three-body terms, in the most general expression for baryon masses. 
We find that powers of M'^ or M'^ do not correlate strictly with powers of the quark 
masses rrid — rriu or — rriu and, as a result, the suggested hierarchy of sizes for terms 
in the most general expression for baryon masses that involve multiple factors of M'^ 
and/or does not hold strictly for dynamical calculations in HBChPT. Numerical 
calculations of the sum rules that give constraints on the three-body terms indicate 

§ The connection of this representation to the usual effective-field methods was discussed in detail in 
[151 [TO] . The results in each case can be summarized in terms of a set of effective interactions that have 
the appearance of interactions between quarks in the familiar semirelativistic or nonrelativistic quark 
models for the baryons [5S1 UDl , but the corresponding matrix elements can be translated back to 
expressions in terms of the relativistic effective baryon fields. 
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that contributions of the three-body hypercharge sphttings are comparable to those 
from the one- and two-body isospin sphttings. In addition, the calculations show that 
contributions of the three-body isospin splitting operators to baryon masses are very 
small. Ignoring terms associated with the three-body isopsin splitting operators, we do 
phenomenological fits to the experimental data using our general expression. Our fits 
provide the values of the parameters that must be explained in the loop expansion in 
ChPT or any other dynamical model. The results of the fits also confirm that three-loop 
contributions will necessarily play a role in explaining masses at the scale of about 1 



Our parametrization has the minimal number of operators needed to describe all 
the octet and decuplet masses and can be translated back into HBChPT [181119]. Our 
general expression for baryon masses is particularly useful to an analysis of the baryon 
mass splittings due to simultaneous hypercharge-breaking and isospin-breaking effects. 

The paper is organized as follows. In Sec. [21 we discuss how the general 
parametrization of the baryon octet and decuplet masses is constructed. A reduction of 
the general expression to the independent one- and two body operators at two-loop level 
and a general proof on the cancellation of the three-body terms are shown in Sec. [3l 
Then, we briefly discuss the hierarchy of sizes of the coefficients of the various operators 
in the general expression for baryon masses, investigate the likely size of the three-body 
terms, and do the phenomenological fits to the experimental data in Sec. HI Concluding 
remarks are presented in Sec. [51 

2. General parametrization of the baryon octet and decuplet masses 
2.1. Construction of the mass operators 

First of all, we want to remind the reader that the operators that appear, while 
nonrelativistic in appearance, are all that are needed to describe the structure of matrix 
elements in exact QCD. 

Before introducing our general parametrization of the baryon masses, we need to 
show how to construct the possible mass operators using the unit operator, 1, and the 
symmetry-breaking factors, M'^ = diag (0,1,0) and = diag (0,0,1), in conjunction 
with the quark spin operators cr. Note that M" = diag (1,0,0) = 1 — M'^ — M'^, so 
symmetry breaking terms in iVf" can be transformed to M'^ and M^. We are only 
interested in the mass operators that can have nonzero matrix elements in or between 
the octet and decuplet states. 

If isospin breaking is ignored, employing the Morpurgo's general parametrization 
method as discussed in [H], one can easily find that the following mass operators 
intervene in the expression of the baryon masses: 



MeV. 
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j^jyfc j^jVfc 

where G u,d,s are flavor indices, 1 and ^j-^^ CTj ■ aj are the flavor-symmetric 

operators, the first line is the result for no flavor symmetry breaking, and the second 
and third lines introduce the hypercharge-breaking operators. Note that the operators 
carry the "quark" labels of the effective fields that indicate the flavor index on which 
they are to act. There are implied unit flavor matrices for the quarks whose labels 
do not appear explicitly. For example, Yli has the complete flavor structure 
{M^nlj/jlk'k + '^i'iMyjlk'k + Ij'jli'i^fc'fc)- III Eq. ([T]), terms in the first and second 
lines are the one- and two-body operators and terms in the third line are the three-body 
operators in the sense that they act nontrivially on one, two, or three indices. 

It is straightforward to generalize the result shown in Eq. ([T]) to find the possible 
isospin splitting operators - terms with one, two, or three factors of M'^. That can be 
done by replacing in turn the matrix by M'^ at every place it appears in a flavor 
splitting operator. For example, the flavor splitting operators M/, ^j-^^ M^Mj, and 
J2i=ij=ik M^MjCTj ■ (Tk generate the foUowings 



YMtM^^^YMf^h E^^^i' Y.MtMf, (2) 

«7^i «7^i «7^i ^T^i 

ii^ji^k i+i+k i+i+k ij^jj^k 

Since the indices i, j, and k of the sums are dummy, one can easily find that some of 
the generated operators are identical. For example, 

Y MfM; = Y M^Mf , Y MtM;cT, ■<T,=Y M^MfcT, ■ a, , 

Y MtM]Ml = Y MtM^Ml = Y MtM^Mt , (3) 

i^j^fc i^jVfc 

and so on. Keeping in mind Eq. ([3]), we find the possible one- and two- body isospin 
splitting operators 

YMf, YMfcT^.a,, YMf^h 

Y MfM;a. ■ '^^^ ' E ' E MfMfa, ■ aj , (4) 

and the possible three-body isospin splitting operators 

Y MfcT, -(Tk, Y MfM;(T, -CTk, Y MfM^a, ■ cr^ , 

iT^jT^k i+j+k i+O+k 

Y MtM]Ml , Y MtM]MlcT, -CT^, Y ^tM]MlcT, ■ cr^ , 
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Y MfMfMt , Y. MfMfMtcT, . a, . (5) 

We note that some isospin splitting operators listed in Eqs. (jl]) and ([5]) (e.g., 
terms involving multiple factors of M*^) are not explicitly shown in the parametrization 
constructed by Dillon and Morpurgo. However, it is worth pointing out that, using the 
identities 

Q = -1 - Af^ - M% g2 = li _ _ (6) 
3 9 3 3 

for the charge operator Q = diag (2/3, — 1/3, — 1/3), and also the matrix relations 

MfMf = MfMf = 0, M'^M'^ = Af^, and M'M' = M% the electromagnetic 

operators introduced by Morpurgo in [15] can be reduced to sums of the flavor-symmetric 

operators, the hypercharge splitting operators shown in Eq. ([T]), and the isospin splitting 

operators (except those with three factors of M'^) listed in Eqs. (jl]) and ([5]), respectively. 

2.2. A general expression of the baryon masses 

Not all the splitting operators shown in Eqs. ([T]), (jl]), and are independent. 
Using the identity Mfi^Mj^ = Mf,-Mj,j, where p = u,d,s, and the exchange operator 
Pij = {1 + (Ti ■ (Tj)/2 from [19] to rearrange indices, we can easily show that 

Mf MJ = MfM^a, ■ cr, . (7) 

As a result, we obtain the following identities 

Y MfM^Ml = Y MfM^Ml^, . a, , (8) 

where p,q = d, s. Taking into account these identities, we find that the numbers of 
operators in Eqs. ([I]), (jlj) and ([5]) are reduced to 8, 5, and 9, respectively. Hereafter, for 
the operators related by Eq. ([8]), we choose to work with those without spin operators. 

We now have 22 operators (2 symmetric, 6 hypercharge splitting, 5 one- and two- 
body isospin splitting, and 9 three-body isospin splitting) to describe a total of 18 
baryon octet and decuplet masses. Therefore, four operators must be redundant. To 
determine what operators are independent, we make the tables of their contributions 
to the baryon masses and then consider the corresponding Gram matrix. For the eight 
flavor-symmetric operators and the hypercharge splitting operators, their contributions 
to the baryon masses are given in Table [1] where the three-body hypercharge splitting 
operators are denoted as 

si= Y ^i"^^ -^k. S2=Y Mi^i^^j -(Tk, ss= Y m:m;m^ . (9) 

i^j^k i+i+k i+i+k 
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Table 1. Contributions of the flavor-symmetric operators and the hypercharge 
splitting operators hsted in Eq. ([1]) to the baryon masses. For simphcity, all the 
sums are suppressed. Contributions of the dependent operators identified by Eq. ([5]) 
are not shown. The matrix elements for all particles in a hypercharge multiplet are 
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We present in Table |2]the contributions of the thirteen isospin sphtting operators where, 
for simphcity, the three-body isospin sphtting operators are labehed as 

ti=Yl MfM-Ml , t2 = 5^ MfM^M^.a, .a,,t,= MfUfa, ■ cr^ , 

t4= Yl MfM^Mt ,t,= J2 MfMfMlai ■a,,t,= Mfufut , (10) 

t7= Y ^i"^^ -(^k, ts= J2 MfM;(T^ -(Tk, tg= Y MfM;(T, ■ CTk • 

iy^jy^k ij^jj^k ii^ii^k 

We then consider the 22 x 22 Gram matrix A^r = M. associated with the 18 x 22 
matrix }A defined by joining the 18 x 8 matrix in Table [T] (extended to all states in 
each hypercharge multiplet) and the 18 x 14 matrix in Table [2l M.^ has a vanishing 
determinant and four zero eigenvalues. This indicates that there are four relations 
among the 22 operators. Applying the conditions for linear dependence of the mass 
operators, A^rX = for the eigenvectors x with zero eigenvalues, we find the four 
following relations: 

= -6 1 + ^ <T, ■ <Tj- + 4 ^ M; - 2 ^ M/cTi ■ CTj- 

+Aj2Mt-2j2Mt(T,-a,-2j2MtM; + 2Y,MfM;cT,-cT,, (11) 

= -6 1 + ^ cr, ■ cr^. + 6 ^ M/ - 2 ^ M/o", ■ cr^- - 2 ^ M/MJ - Si + 2s2 

+4 ^ - 2 ^ Mf(Ti ■ cr, - 4 ^ MfM] + 2ts + 2tg , (12) 

= -6 1 + ^ <Ti ■ o-j + 6 ^ M/ - 2 ^ M/o-i ■ o-j - 2 ^ M/M| - si + 2s2 
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Table 2. Contributions of the isospin splitting operators listed in Eqs. (jl]) and ^ to 
the baryon masses. Here, for simpHcity, we denote ci — J2i ^^ti ^2 — J2i^j ^iff^i ' ""j) 
C3 = T.^^J MfMl C4 = E.^j MfM^a, ■ a^, and cg = Y.^^^ MfUf. Contributions of 
the dependent operators identified by Eq. ([H]) are not shown. 
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+2 ^ M/ - 2 ^ MfcTi ■ o-j - 2 ^ M; + 2 ^ Mf 

+t7 + 2^8 - 2t3 - 2^4 + 2^5 , (13) 

= 2j]M/ + 2^ M/M; + Si - 2s2 + 2 ^ M/m; - 2^8 - 2ti + 2t2 . (14) 

Using these relations, we can select 4 operators to be dependent. The chosen 
operators are J^i^j ^f^j^i ' ^jy tr^ ts, and tg. Note that the chosen two-body 
operator is related in Eq. f|TT]) only to other one- and two-body operators and does 
not change the values of any three-body operators. On the other hand, it follows from 
the remaining equations that the matrix elements of ty, ^8^ and tg can be absorbed 
in the matrix elements of one- and two-body operators and the remaining three-body 
operators (si — 2s2), ti to t^. No one- or two-body operators are mixed into the three- 
body operators, so the expected smallness of the remaining three-body terms should be 
maintained. 

We now have a following set of 18 independent mass operators 
Flavor-symmetric 1 , cTj ■ ctj , (15) 

Hypercharge splitting 

One and two-body: ^ , ^ M/o-^ ■ aj , ^ M/M; , 
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Three-body: J^M^a.-ak, M^M^a, ■ ct^ , J^M^MJM',. (16) 

iy^jy^k ii^ii^k i+i+k 

Isospin splitting 



One and two-body: ^ M"' , ^ Mf o-^ ■ o-^- , ^ MfU. 



d 
j 



Three-body: ^ M/M/o-^ ■ cr^ , ^MfMfMt, (17) 

Mixed sphtting 

Two-body: ^ , 

Three-body: ^ M/m;M,V, ■ o", , J^MfMfM',, 

MfMfMlcT^.cT,. (18) 

Note that the mixed sphtting operators affect both hypercharge and isospin sphttings 
and involve at least one factor each of and M'^. 

We can now construct a general expression for baryon masses. Using the flavor- 
symmetric operators and other independent mass operators given in Eqs. (fT5|) - (fT8|) . 
we write the most general expression for baryon masses as 

'Hb = + A ^ (Ti- (Tj + Hlnter + 'Hintra , (19) 

where TCinter and H-intra are the general expressions for the intermultiplet splittings and 
intramultiplet splittings, respectively. Namely, 

Hlnter = 5 ^ M/ + C ^ M/cT, ■ CT, + D ^ MfM/ + aSi + hS2 + CS3 , (20) 

and 

6 

Hlntra = A' Mf + B' ^fcT, ■ CT, + C ^ + D' ^ + "^^^^ ' ^^l) 

i iytj iytj iytj 1 = 1 

Here mo. A, B, C, D, A', B', C, D', a, b, c, and rfj {i = 1, ...,6) are the parameters. 



3. A parametrization of baryon masses to two loops 
3.1. Generation of the one- and two-body operators 

In [19], we analyzed the structure of meson loop corrections to the 0{ms) expressions 
for the baryon masses. Our approach was based on the standard chiral Lagrangian of 
the heavy-baryon chiral effective field theory [21]. However, we used a spin- and flavor- 
index or "quark" representation of the effective octet and decuplet baryon fields rather 
than the usual matrix expressions for the fields. There, to one loop for equal mass light 
quarks in an expansion in m^, we found the appearance of all one- and two-body mass 
operators listed in Eq. ([T]), but no three-body operators. 
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For the case of nonzero — m„, it is straightforward to show that one-loop 
corrections involving M'^ insertions and the effects of meson mass differences in loop 
integrals introduce all possible one- and two-body isospin splitting operators shown in 
Eq. (jl]). We find no appearance of the three-body isospin splitting operators at this 
level. 

Recall that, in [20], our calculations of the electromagnetic contributions to 
baryon masses were carried out including the one-loop mesonic corrections to the basic 
electromagnetic interactions, so to two loops overall. To this order, the electromagnetic 
contributions also produce a complete set of the one- and two-body mass splitting 
operators, but again no three-body operators are generated. A general proof on the 
cancellation of the three-body terms up to two loops will be given in the next subsection. 

As a result, a general parametrization of baryon masses to two loops is 

nB = mol + Aj2(^i-(^j + Bj2M'+Cj2 ^t^^ -^i + dY^ M^M; + n'lB , (22) 

where Ti'iB is the isospin splitting mass Hamiltonian consisting of the one- and two-body 
operators only, 

n'lB = A'Y,Mf + B'Y, Mf(Ti ■ (Tj + C'Y, MfM; + D'Y, Mfuf . (23) 

Again, mo. A, B, C, D, A', B', C, and D' are the parameters. 

Before moving on to the next subsection, we want to make clear that, in this 
subsection, we are not deriving the parametrization from the calculations, but are 
summarizing what happens in actual dynamical calculations. We emphasize that the 
form of the parametrization above follows rigorously from properties of QCD as shown 
by Morpurgo [11], and is encompassed also in the most general chiral Lagrangian in 
effective field theory. While we obtained the same one- and two-body operators in our 
earlier calculations in HBChPT [181 [13 EQ] as summarized here, the form of the general 
parametrization, whether approached from QCD or effective field theory, is independent 
of the results of particular dynamical calculations. The situation is different in the next 
subsection, which deals specifically with the the origin of three-body effects in HBChPT. 

3.2. Cancellation of the three-body terms 

We present here a general proof on the cancellation of the three-body terms through two 
loops in HBChPT. Our analysis is done using an expansion in meson (or photon) loops 
in time-ordered perturbation theory. The four possible types of time-ordered diagram 
that involve all three indices i, j, k ("three-body operators") and one- or two meson (or 
photon) loops are shown in Fig. [H We are going to show that these diagrams cancel 
exactly with the renormalization diagrams with the same topology. 

In Fig. [H a solid vertical line represents a "quark" (a flavor index) moving upward 
in time from an initial state with flavor index i to a final state with index i' rather than 
a propagating QCD quark. These flavor indices i,j, k and k' are connected in the 
effective field theory to the corresponding indices on initial and final effective baryon 



A parametrization of the haryon octet and decuplet masses 



11 



fields ip^ji^, '4'i/' j' k' described in detail in [THl [12], and the spin matrix elements are 
calculated with respect to these fields. The baryon fields are suppressed in the present 
graphical notation, and only the connections of initial and final fiavor indices and the 
structure of the spin-dependent operators are shown explicitly. 

A bold solid line connecting fiavor lines represents the exchange of a particle in the 
chiral effective theory such as a meson or photon between fiavor (or "quark") lines. This 
exchange corresponds in the baryon picture to a diagram with a meson or photon loop 
connecting the baryon to itself, with the vertex operators acting on different flavor and 
spin indices at different times in time-dependent perturbation theory. The counting of 
loops is done in this baryonic context. Note that the horizontal dot-dashed lines pick 
out the intermediate states. The vertex operators typically involve flavor operators such 
as Gell-Mann A matrices or charge or mass matrices, and spin operators expressed in 
terms of Pauli matrices. 

Similarly, a squared dot represents an operator or mass insertion that acts only on 
one index, while a zigzag line represents a "two-particle" or two-index insertion with no 
propagator or time ordering. 

Note that in Figs. [H^a), [11(c), and[I]^d), only one of the possible time orderings 
is shown and in Fig. [11(b) both the possible orderings are included. There are two, 
sixteen, and again sixteen possible orderings for each type of diagram in (a), (c), and 
(d), respectively. 

Firstly, let us consider Fig. [Tl(a) depicting a one-loop exchange diagram with a one- 
particle "mass" insertion. As will be illustrated below, the value of the diagram can be 
expressed as a product of energy denominators and a matrix element of vertex operators 
dependent on the flavor and spin indices. In the flgure, the particle exchange connects 
the quark lines j and k and the "mass" insertion is on the quark line z at a time between 
two vertices of the exchange line. The "mass" insertion could include a mass matrix, 
a spin, a charge matrix, a moment operator, and so on; the latter two appeared in our 
analysis of baryon moments. The matrix element involves all three indices i, j, k so the 
diagram is three-body in Morpurgo's classiflcation. However, the three vertices are on 
separate lines so they are not constrained topologically. Moving the j and k vertices 
along those quark lines to above or below the i vertex, which can be done freely, gives 
diagrams topologically equivalent to renormalization diagrams, with the same vertex 
factors. In addition, the energy denominator for the diagram corresponds to the one 
in the renormalization constant associated with the exchange. Multiplying the single 
particle mass insertion amplitude, which corresponds diagrammatically to Fig. [H^a) 
with the particle exchange deleted, by the renormalization constants for the initial and 
flnal states associated with the particle exchange on the other lines gives a term with the 
same energy denominator. This cancels the three-body contribution from the diagram 
under consideration. 

As a demonstration, we now show the time-ordered perturbation theory calculations 
in HBChPT for the diagram in Fig. [H^a) when the exchange is a meson. This 
corresponds in the baryon picture to a one-loop diagram with the meson loop connecting 
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Figure 1. Time-ordered one- and two- loop three-body diagrams which cancel exactly 
with renormalization diagrams. A solid vertical line represents a quark moving upwards 
toward later times. The horizontal dot-dashed lines pick out the intermediate states. 
A squared dot indicates a one-particle or one-index insertion. The zigzag line is a two- 
particle or two-index insertion with no propagator or time ordering. A bold solid line 
represents the exchange of a particle in the HBChPT Lagrangian between the baryon 
and itself (a loop diagram in the baryon picture), (a): One-loop exchange diagram 
with an one-particle "mass" insertion, (b): The exchange diagrams with a two-particle 
"mass" insertion where both time orderings of the particle exchange vertices are shown, 
(c): A double exchange diagram with the exchanges connecting the quark lines i and 
and j and k. (d): A double exchange diagram similar to (c), but with an one-particle 
insertion on one of the quark lines. Note that in (a), (c), and (d), only one of the 
possible time orderings is shown and in (b) both the possible orderings are included. 
There are two, sixteen, and sixteen possible orderings for each type of diagram in (a), 
(c), and (d), respectively. 
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the baryon to itself. An energy denominator 1/{Eq — En) appears in the expression for 
the perturbed energy for each intermediate state \n) of the baryon system between 
successive vertices. In the heavy-baryon approximation, any internal baryon momenta 
k resulting from baryon interactions are small on the scale of the baryon mass M and 
can be neglected [21]. The baryon energy is then always M for the initial baryon 
at rest, and cancels out in the difference Eq — En- Hence, the energy factor reduces 
simply to —^/YliJ^i^ where the sum is over the energies of the lines cut in the 
intermediate state, but with no contribution from the quark lines. In this example, 
E, = Ei{k) = y/fc2 + Mf, where Mi is the mass of the meson in the loop and k is its 
momentum in the baryon rest frame. We find that a contribution from the diagram 
shown in Fig. [T](a) is 

Tpl^ / EAkW ^'' i>^f,^rkcT,-kXl,) 



(27r)3 2y'k^ + A.f^ 

= i5^/;a,.(A;,,.^,-^fcAy , (24) 

I 

where Oj/j is the "mass" insertion vertex, A is a Gell-Mann matrix in flavor space, E{k'^) 
is a form factor used to regularize the integral, and // is a modified integral defined as 

A contribution from the renormalization diagram can be found by multiplying 
the "mass" -insertion amplitude Oi/i by the wave function renormalization constants 
Z = 1 — 6Z for the initial and final states associated with the meson exchange on line 
j and k. The result is 

— 2 ^^fk' jk + ^Zjiy.jk Oi'i) = —Oi'i dZyy.jk , (26) 

where 

SZj>k>.jk = 2 ^^j'j^k'k (^j (^k- (27) 

I 

It is obvious that contributions from the diagram in Fig. [T]^a) and from its 
renormalization diagram cancel exactly as suggested by its topology. 

Secondly, we turn to the case of the two-particle "mass" insertion. Again, we 
consider the exchange diagrams with a particle exchange connecting the quark lines j 
and k, but with a two-particle "mass" insertion "connecting" quark lines i and j as 
shown in Fig. [It^b). Note that the two-particle insertion, denoted by Oi>ij>j, is an 
instantaneous short- distance term (e.g., OiUjij = {M^cr ■ cr)j/jj/j, where p = d,s) but 
with one-loop structure, so we count it as equivalent to one loop. For the first (second) 
ordering in Fig. [U^b), since Oi'ij'j has no effect on the vertex on line k, and its vertex 
on line j is above (below) the exchange vertex on that line, the exchange vertices can 
be moved freely down (up) on lines j and k to obtain a figure with the topology of a 
renormalization diagram. On the other hand, because the two-particle insertion has no 
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energy denominator and does not contribute to 6Z, the energy denominators for both 
the orderings in Fig. [I](b) are identical to those for the renormahzation diagrams with 
the same topology. As a result, when all contributions are added together, the three- 
body terms from the exchange diagrams and their renormahzation diagrams cancel each 
other as expected. For example, if the exchange particles are the mesons, the diagrams 
in Fig. [It^b) give the following contribution 

o^^'l {^i'id'j"^r'j^j'^k>^i'k + ^fj"(^j"-{(^i'i,j"j(^k>^i'k)) ) (28) 

that cancels exactly with the contribution from the renormahzation diagrams 

~ {C^i'iJ'j" ^Zj/ik'-Jk + SZj'k';j"k ; (29) 

where SZ is given by Eq. (!27|) . 

Thirdly, in Fig. [I](c), we consider a double exchange diagram with the exchanges 
connecting, say, ij and jk. The vertices on lines i and k can be freely slid along those 
lines, while the two vertices on j can be slid in opposite directions. Moving the ij 
exchange up leads to a figure topologically equivalent to that encountered for an ij 
exchange multiplied by the renormahzation constant for a jk exchange, or conversely. 
Recall that there are sixteen possible time orderings for this type of diagram. When we 
combine the energy denominators for various orderings, we in fact get those for products 
of single exchange diagrams and renormahzation factors. We also need to include 
explicitly the contributions of single ij and jk exchanges multiplied by renormahzation 
factors. The cancellation occurs only when all contributions are added together. 

As an example, we take the exchanges to be the mesons and find the following 
contribution from all the time orderings of the double exchange diagram 
2 

- - ^ {Iil'i' + I'lh') Vi'j'k'-ijk , (30) 
1,1' 

where Ii is the common integral for the exchange of meson /, 

f}2 poo 7,4 

Ii= / dk—^ oF\k''), (31) 

the modified integral I[ is given by Eq. (!25l) . and the group factor Vi'j'k'-ijk is 

Vi'j'k';ijk = (cTi ■ <^j){<^j ■ <^k) (^A[,j(A'A' )j'j\^k'k + K'ii^^ ^^)i'j^i'k) ■ (32) 

To determine a contribution from the renormahzation diagrams, the amplitudes for the 
two separate ij and jk exchanges are multiplied with the parts of the renormahzation 
constant Z associated with other exchanges. Sum of the two products gives terms that 
cancel the original "three-body" amplitude shown in Eq. fl5Ul) . 

Fourthly, it is straightforward to generalize the above arguments for a double 
exchange diagram with the exchanges connecting two pairs of lines and an extra one- 
particle "mass" insertion as depicted in Fig. [T](d). Since we can always slide one 
exchange line topologically outside the rest of the diagram, whatever the initial time 
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Figure 2. Three distinct topologies of the three- loop three-body diagrams that cannot 
be cancelled by renormalizations. 





orderings of the vertices, the cancellation occurs when all time orders are summed and 
all the (5Z's associated with that line in the various time-orders are included. 

Finally, we present in Fig. [2] the three-exchange diagrams that cannot be cancelled 
by renormalizations. The same three cases also appear when one of the exchanged 
lines is replaced by an instantaneous two-particle "mass" insertion. One can see that 
there are three distinct topologies involved: i) two vertices on each line with the loops 
entangled; ii) one vertex on each of two lines and four vertices on the remaining line 
with a loop from that line to itself enclosing the other two vertices; and iii) two vertices 
on each of two lines with crossed exchanges, and a third trapped vertex on one of 
the lines with the exchange connecting to the third line. In the remaining three-loop 
diagrams, one line can always be slid outside the remainder of the diagram, leading 
topologically to a renormalization-type diagram, and the final result is actually cancelled 
by renormalization terms. 

We conclude this subsection by noting that since the three-body contributions only 
appear at three loops, one can expect their contributions to the baryon masses to be 
very small. 



4. Hierarchy and Fits 
4.I. Hierarchy of sizes 

In the case of hypercharge breaking, Morpurgo [23l [2l] argued on the basis of ideas on 
the interactions of gluons with structure (rather than elementary) quarks and fits to 
some of the coefficients that the sizes of the coefficients of the various operators shown 
in the most general expression for baryon masses (see Eqs. f|T9|) - fl2Tl) ) should satisfy 
a set of hierarchical relations, with suppressions by a factor of ~ 1/3 for each factor of 
M'^ and by a factor in the range from 0.22 — 0.37 for each extra gluon exchange needed 
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to get the number of indices and spin factors in the operator [TBI EB] (l- 

Somewhat different arguments apply in the HBChPT approach when the higher- 
order operators are generated by meson loop corrections to a basic set, as in [21] and |19j . 
Factors of then appear multiplying differences of loop integrals involving kaons or 
eta mesons and pions. The coefficients are small because these differences would vanish 
for equal-mass mesons. also appears in explicit non-calculable mass insertions with 
fairly small coefficients, smaller for spin-dependent than spin-independent insertions. 
These effects were studied in detail in our previous work |T9] . 

Similar results hold in the case of isospin breaking, with factors of M'^ now 
multiplying either differences of pion (or kaon) loop integrals for mesons within the 
isospin multiplet, or explicit c?-quark mass insertions. Since << rris, and the mass 
differences within the pion and kaon multiplets are much smaller than those between 
different multiplets, the coefficients of M'^ are expected to be much smaller than the 
corresponding coefficients of M*. The coefficients are also expected to be smaller for 
spin-dependent than for spin-independent forms. The results are consistent with the 
observed accuracy of various sum rules for the masses. 

It is important to note that the sizes of terms proportional to M'^M'^, M'^M'^, or 
M^M^ do not go strictly as (m^ — my)"^, {rria — mu){ms — rriu), or {nis — m^)'^ as mass 
insertion arguments would seem to suggest. The first two can both be proportional 
instead in the loop expansion of HBChPT to rrid — rriu because of the structure of the 
loop corrections, while the last can vary as — m„. This can be seen explicitly in the 
M*cr ■ cr and M'^M^cr ■ cr terms obtained from the one-loop mesonic corrections to the 
baryon masses in [T9|. The leading M*-dependent corrections, given in Eq. (3.14) in that 
paper, are all proportional to differences of meson loop integrals which vanish as m^— 
as the quark and meson masses become equal. A single factor of M would therefore be 
expected. However, because of the particular tensor structures involved in the products 
of Gell-Mann A matrices in the two-body couplings, the loop corrections that involve 
meson exchanges between quarks contribute to both the M^cr ■ cr and M^M^cr ■ cr terms, 
with similar sized corrections. The first is similar in size to other contributions to the 
M*cr ■ cr coefficient. The second is much larger than would be expected from the simple 
QCD argument. 

Similar results hold for the meson loop corrections that involve isospin breaking. 
Since terms in M'^M'^ and M'^M^ can be of the same order as terms with a single factor 
of iVf^, it is not justified in applications to drop them without further analysis. Thus, 
while contributions that explicitly involve higher powers of the quark or meson mass 
differences, or higher loop integrals, are clearly (or presumably) suppressed, we do not 
expect to find Morpurgo's strict hierarchy for terms that involve multiple factors of M'^ 
and M^ 

II Dillon and Morpurgo have adopted the value of 0.3 for the suppression factor due to each extra gluon 
exchange [28] . 
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4- 2. Fits to the data 

Before doing the fits, we want to study the constraints on the three-body terms. Note 
that without the three-body terms, we have nine parameters mo. A, B, C, D, A', B', 
C, and D' to describe eighteen mass states of the baryon octet and decuplet. Hence, 
we expect there to be nine sum rules among the baryon masses. To find these sum 
rules, we construct a 18 x 18 7V1 matrix by replacing all columns associated with the 
three-body operators in the 18 x 18 7V1 matrix defined earlier in Sect. 12.21 with zero 
entries. The obtained matrix has nine zero eigenvalues. The sum rules are just the 
inner products of the left null eigenvectors Xq with /i = Xiv, XqH = Xq/Av = 0, where 
f is a eighteen-component column vector of coefficients mo. A, B, C, D, a, b, c. A', B', 
C, D', and di (i = 1, 6). We find 

A° - A+ = n - p 

A- - A++ = 3{n-p) 

A° - A++ = 2{n -p) + (S° - S+) - (S- - S°) 

5- -S° = (S--S+)-(n-p) (33) 
E*- - S*° = (S*- - -{n-p) 

2S*° - S*+ - S*" = 2S° - S+ - 

Q- - A++ = 3(S*° - 

(3p - n)/2 + S° - (3A + S+ + S° - ^-)/2 = (3H*° - S*~)/2 - 2S*+ + A++ . 

The first six relations are the well-known sum rules for isospin splittings. The fourth is 
the Coleman-Glashow relation, suggested originally on the basis of an unbroken SU(3) 
flavor symmetry [l] (another original SU(3) sum rule is the Gell-Mann - Okubo (GMO) 
mass formula [25l|26]). All the sum rules were later established for nonrelativistic quark 
models with only one- and two-body interactions independently of the flavor symmetry 
breaking [21 El HI [5] . Note that the ninth relation is the modified GMO sum rule that 
first appeared in the quark model context [3]. If there are no intramultiplet splittings, 
it reduces to the usual GMO rule. 

The above sum rules are violated by three-body terms. If we transfer all the terms 
to the left hand sides of the equations and denote the differences by 5j (i = 1, 9), the 
results for these sum rules in terms of the three-body coefficients are 

6i = 6(^3 , 62 = 18(^3 + 6(^6 , S3 = Gds + 2(^4 — Ad^ , 

64 = 2di - Ad2 - 24 - 4^5 , 65 = 2di + 2^2 - 64 - 24 - 24 , = -64 - 6^5 , (34) 

6t = — 12a — 6b , Ss = 6c, 69 = —66 + rfi + 4 + 4 — 24 • 

Isospin splittings: We have mentioned at the beginning of this section that the 
coefficients of the mass operators in Eq. f|T9l) are expected to satisfy an approximate 
hierarchy of sizes. As discussed below, the sizes of some three-body coefficients can be 
estimated by evaluating the 5's using the experimental values of the accurately known 
baryon masses. 
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The A baryon masses are not determined with sufficient accuracy for the first three 
sum rules to give a real test of this expectation. The results for the next three as written 
are, in order, ^4 = -0.31±0.25 MeV, ^5 = 0.09±0.93 MeV, and Sq = -1.06±1.18 MeV, 
all consistent with zero within the experimental uncertainties. No significant violations 
of the sum rules are evident. 

Note that the coefficients di {i = 1,...,6) associated with the three-body isospin 
splitting operators are expected to be very small. In particular, the coefficient d^ of the 
three-body operator t^ = J2ij^jj^k M^M^Mj^ that contributes only to the A~ mass, is 
expected to be very small since carries three factors of M'^ (its size is proportional 
at least to (m^ — m^)^). We also expect that the coefficients d^ associated with the 
three-body operators ^5 = ^ ..^^y^, MfMjM^CTi ■ cr^ should be smaller than the other di 
{i = 1, 4) since t^ involves two factors of M'^, one factor of M*, and a spin interaction 
with the third particle. Therefore, we will ignore d^ and d^ when trying to estimate the 
other coefficients d^. Using the values of ^4, ^5, and 6q, we find di — d^ = 0.33±0.50 MeV, 
d2 = 0.24 ±0.25 MeV, and 4 = 0.18 ±0.20 MeV. These results show that contributions 
of the three-body isospin splitting operators to baryon masses are indeed very small and 
can be ignored. 

Hypercharge splittings: If we set di = {i = 1, ...,6), as the above results 
suggest, we are left with the relations 

67 = -12a -6b, §8 = 60, Sg = -66 . (35) 

To estimate the coefficients of the three-body hypercharge splitting terms a, b, and c, we 
use the average mass of the A's and other accurately known baryon masses to evaluate 
6y, 6s — Sg (note that A^+ does not appear in this combination), and ^g. The result 
are 67 = -23.54 ± 0.55 MeV, §8-5^ = 3.40 ± 1.08 MeV, and = -9.95 ± 2.24 MeV. 
Then, Eq. ([35]) gives a = 1.13 ± 0.06 MeV, b = 1.66 ± 0.37 MeV, and c = -1.09 ± 0.41 
MeV. These contributions of the three-body hypercharge splitting terms are statistically 
significant and as large as observed isospin splittings, so they should be kept in a 
complete analysis of the baryon octet and decuplet masses. 

General constrained fit: We now consider a ffi to the baryon masses using the 
general expression in Eq. (fT9l) neglecting the coefficients di {i = 1,...,6). Using 12 
parameters mg. A, B, C, D, A', B', C, D', a, b, and c, we can do a weighted least- 
squares ffi to the 15 measured quantities (14 accurately known baryon masses plus the 
average mass of the A's). 

A best ffi is obtained at values (in MeV) 

mo = 1082.9 ± 1.0, 

A = 24.66 ± 0.17 , B = 183.66 ± 1.05 , C = -17.10 ± 0.34 , D = -2.94 ± 0.71 , 
a = 1.22 ± 0.17 , 6 = 1.52 ± 0.35 , c = -0.93 ± 0.40 , 

A' = 0.95 ± 0.18 , B' = -0.60 ± 0.08 , C = 1.69 ± 0.24 , D' = 0.78 ± 0.04 . (36) 

The ffited masses have an average deviation from experiment of only 0.07 MeV and a 
= 1.73 (with 3 degrees of freedom). The calculated values of the baryon masses are 
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Table 3. Baryon masses in units of MeV. Excluding the A resonances, the average 
deviation |AMs| = 0.07 MeV, where AMb = Mg^^^^^ - M^^P*'. The experimental 
data are from the listings given by the Particle Data Group [27]. No definite masses 
or uncertainties are given for the A resonances as different experimental analyses are 
in conflict. Only ranges of mass are given here. 



Baryon 


Theory 


Expt. 


lAMisI 


P 


938.27 ± 1.50 


938.27 ± 0.00 


0.00 


n 


939.57 ± 1.51 


939.57 ± 0.00 


0.00 


A 


1115.68 ± 2.10 


1115.68 ± 0.01 


0.00 


E+ 


1189.39 ± 2.28 


1189.37 ± 0.06 


0.02 


1]° 


1192.64 ± 2.30 


1192.64 ± 0.02 


0.00 




1197.45 ± 2.36 


1197.45 ± 0.03 


0.00 




1314.64 ± 3.58 


1314.83 ± 0.20 


0.19 




1321.39 ± 3.62 


1321.31 ± 0.13 


0.08 


A 


1232.00 ± 1.48 


1232.00 ± 2.00 


0.00 


A++ 


1230.82 ± 1.46 


1230.5 ± 0.2 - 1231.88 ± 0.29 




A+ 


1230.57 ± 1.47 


1231.6 - 1234.9 ± 1.4 




AO 


1231.87 ± 1.53 


1233.1 ± 0.3 - 1234.35 ± 0.75 




A- 


1234.73 ± 1.64 








1382.74 ± 1.95 


1382.80 ± 0.40 


0.06 




1384.18 ± 2.02 


1383.70 ± 1.00 


0.48 




1387.18 ± 2.06 


1387.20 ± 0.50 


0.02 




1531.81 ± 3.37 


1531.80 ± 0.32 


0.01 




1534.95 ± 3.41 


1535.00 ± 0.60 


0.05 


n- 


1672.45 ± 6.78 


1672.45 ± 0.29 


0.00 



given in Table [3] where the experimental data are from the listings given by the Particle 
Data Group [27J. No definite masses or uncertainties are given for the A resonances as 
different experimental analyses are in confiict. Instead, we give here their mass ranges 
except for A~. Our predictions for the A resonances agree with the experimental values 
within the margin of error. 

Note that the best-fit values of the parameters A', B', C, and D' show that terms 
in M'^M'^ and M'^M^ are of the same order as terms with a single factor of M*^, as 
expected in the meson loop expansion, rather than much smaller as suggested by the 
QCD hierarchy arguments. 

Baryon mass differences: We remark that one can use the general expression for 
baryon masses in Eq. (fT9l) to describe the baryon mass differences. Then the parameters 
associated with two flavor-symmetric operators and with the intermultiplet splitting 
operators drop out, so there remain 10 parameters {A', B', C, D', and di {i = 1, 6)). 
Using contributions of the isospin splitting operators to the mass differences given in 
Table [21 we easily find the foUowings 

n-p = A' + 2B' + 2D' - 4^ , 

S- - S+ = 2{A' ~ B' + C' + D' - 2d3 + 4 - 2d5) , 

E" - E° = A' -B' + C' + 2D' - Ads + 2^4 - 4^5 , 
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A' - AB' + 2C' + 2di - Ad2 , 

2{A' + 2B' + C' + D' + d3 + d4 + 4) , (37) 
A' + 2B' + C' + 2D' + 24 + 24 + 24 , 
A + 2B' + 2C' + 2di + 24 , 
2{A' + 2B' + D' + 4) , 
3(A' + 2B' + 2D' + 24 + 24) , 
A' + 2B' + 2D' + 24 . 

If the three-body isospin sphtting terms 4 are neglected, we have 4 parameters 
A\ B', C, D' to describe the baryon mass differences. We studied this case earher in 
pm [22] using the one- and two-body electromagnetic operators. A best weighted fit to 
the seven known mass splittings other than those for the A baryons is obtained at values 
(in MeV) oi A' = 0.91 ±0.19, B' = -0.59 ±0.08, C" = 1.73 ±0.23, and D' = 0.78 ±0.05 
with an average deviation from experiment of 0.13 MeV and a = 1-67 (with 3 degrees 
of freedom) . The calculated values of the baryon mass splittings are the same as those 
shown in Table I in |22j . 

Our above fits provide the parameters that must be explained in the loop expansion 
in chiral effective field theory or any other dynamical model, and show that three-body 
contributions will necessarily play a role in explaining masses at the scale of about 1 
MeV. The scale for three-body isospin effects is at least ten times smaller. 

5. Conclusions 

We have constructed the general parametrization of the baryon octet and decuplet 
mass operators including three-body terms using the unit operator and the symmetry- 
breaking factors, Af^ and M*, in conjunction with the spin operators. Our general 
expression for baryon masses includes all order of the light quark mass difference 
"^d — ^u- Using the Gram matrix analysis, we established a minimal set of 18 
independent operators to describe 18 baryon octet and decuplet mass states. 

We have demonstrated the cancellation of the three-body terms through two loops 
in HBChPT and have identified three distinct topologies of the three-loop diagrams 
that generate the three-body terms. We also investigated the likely size of the three- 
body terms using the sum rules that give constraints on their coefficients. Numerical 
calculations of the sum rules indicate that contributions of the three-body hypercharge 
splittings, while small, are comparable to those from the one- and two-body isospin 
splittings. The contributions of three-body isospin splitting operators are statistically 
consistent with zero, and these operators can be neglected at the present level of 
experimental accuracy. We find that the suggested hierarchy of sizes for terms that 
involve multiple factors of Af^ and/or may not be evident because powers of M'^ 
and do not correlate strictly with powers of the quark masses ma — mu and ms — mu, 
respectively. 



:*0 



AO-A+ 
A- - A^ 
AO-A+ 
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We have done the phenomenological fits to the experimental data using our general 
expression, but ignoring terms associated with the three-body isopsin splitting operators. 
Our fits provide the values of the parameters that must be explained in the loop 
expansion in ChPT or any other dynamical model. The results of the fits also confirm 
that three-loop contributions will necessarily play a role in explaining masses at the 
scale of about 1 MeV. 

Finally, our parametrization has the minimal number of operators needed to 
describe all the octet and decuplet masses. It can be translated back into the language 
of the heavy-baryon effective field theory and is particularly useful to an analysis of the 
baryon mass splittings due to both hypercharge-breaking and isospin-breaking effects. 
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